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ABSTRACT 

The stellar initial mass function (IMF) is a key property of stellar populations. There is growing evidence that 
the classical star-formation mechanism by the direct cloud fragmentation process has difficulties reproducing 
the observed abundance and binary properties of brown dwarfs and very-low-mass stars. In particular, recent 
analytical derivations of the stellar IMF exhibit a deficit of brown dwarfs compared to observational data. Here 
we derive the residual mass function of brown dwarfs as an empirical measure of the brown dwarf deficiency 
in recent star-formation models with respect to observations and show that it is compatible with the substellar 
part of the Thies-Kroupa IMF and the mass function obtained by numerical simulations. We conclude that the 
existing models may be further improved by including a substellar correction term that accounts for additional 
formation channels like disk or filament fragmentation. The term “peripheral fragmentation” is introduced here 
for such additional formation channels. In addition, we present an updated analytical model of stellar and sub¬ 
stellar binarity. The resulting binary fraction and the dynamically evolved companion mass-ratio distribution 
are in good agreement with observational data on stellar and very-low-mass binaries in the Galactic field, in 
clusters, and in dynamically unprocessed groups of stars if all stars form as binaries with stellar companions. 
Cautionary notes are given on the proper analysis of mass functions and the companion mass-ratio distribution 
and the interpretation of the results. The existence of accretion disks around young brown dwarfs does not 
imply that these form just like stars in direct fragmentation. 

Subject headings: binaries: general — brown dwarfs — methods: numerical — methods: statistical — stars: 
low-mass — stars: luminosity function, mass function — 


1. INTRODUCTION 

The stellar initial mass function (IMF) is a key tool for 
star-formation research because it mirrors the processes of the 
formation of stellar populations (Bastian et al.|2010||Kroupa 
et al. 2013|l. Consequently, the IMF has been subject to exten¬ 


sive research both observationally and theoretically. In recent 
years the majority of the star-formation community has fa¬ 
vored the assumption of continuous star formation from the 
lowest-mass brown dwarfs (BDs) to the most massive stars 
(]Padoan & Nordlundl|2002| |2004| |Hennebelle & Chabrier) 
2008)1. However, there is evidence for at least two separate 


formation channels for most BDs on the one hand and most 
stars on the other. For instance, careful analysis of the ob¬ 
servational data reveals a disagreement between the theoret¬ 
ical predictions of the binary separation distributions of BDs 


and stars and the observed one s ( Bouy et al.|2003| Burgasser 
|et al.|2003||Martm et al.|2003| ana Tclose et al.|2003| ). In this 
light the key assumption of a uniform or continuous formation 
mechanism assumed in most star-formation models needs to 
be questio ned. Based on o bservational data and N-body com¬ 
putations, |Parker & Goodwin| (|2011 j) conclude that the birth 
population of very-low-mass binaries with system masses be¬ 
low 0.2 Mq must be very different from that of M-dwarfs. 


Another important issue is the “BD desert” (McCarthy, Zuck- 
erman, & Becklin |2003 i, an observed dearth of B D compan¬ 
ions to stars. Unli ke earlier studies like e.g. by Gret her &] 
|Lineweaver| (J2006) who interpret this as being related to the 
compamon mass ratio rather than the ab solute mass, the more 
recent survey by Dieterich et al.](]2012j! deduce a lower abso¬ 
lute mass limit o 


companions to stars close to 0.1 Mq. 
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Analytical approaches by jPadoan & Nordlund| ( 2002 
PN02, from here) and [Hennebelle & Chabrier| ( ]2(") 08 2009 


HC08, and HC09, hereafter) deduces the IMF from an analy¬ 
tical description of the distribution of prestellar cloud cores. 
While the stellar IMF is reproduced by these approaches, a 
significant deficit of BDs and very-low-mass stars (VLMSs, 
between about 0.08 and 0.2 Mq) with respect to observa¬ 
tionally constrained IMF descriptions ( [Kroupa|2001 Chabrier] 
2005 1 Thies & Kroupa| [2007 2008] the latter three from 
here on referre d to as C 05, TK07, and TK08, respectively, 
and [Kroupa et al.|2013) appears if realistic properties of the 
prestellar clouds"are assumed. The fundamental reason why 
direct fragmentation of a turbulent molecular cloud rarely pro¬ 
duces BDs is that the formation of a BD requires a high- 
density gravitationally self-bound but very low mass fluctu¬ 
ation that cannot draw significant amounts of additional mass 
from an accretion reservoir (see also |Adams & Fatuzzo|1996| ). 
However, HC08 and HC09 speculate that this deficit might 
be solved by a refinement of their models by including tur¬ 
bulent fragmentation descriptions. A different interpretation 
is that their model is overall correct but reveals that an ad¬ 
ditional formation channel is required to match the observa¬ 
tions. The IMF by C05 is an empirical and almost equivalent 
update of the IMF by |Chabrier| (|2003|>. Fully hydrodynam- 
ical comput atio ns of whole s tar-forming clouds by e.g. |Bate| 
et al. ( 2003] ) and Bate ( 2009aj ) reproduce the formation of BDs 


largely from fragmenting circumstellar disks, although these 
simulations tend to overproduce BDs unless radiative heating 
is included ( |Bate|2009b |2012] >. 

In this paper we introduce the residual mass function 
(RMF) as a correction term for this BD deficiency in the ana¬ 
lytical approaches by PN02 and HC08 with respect to the 
observation-based IMFs in C05 and TK07. In Section |3 the 
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model is described and the RMF is defined. The results are 
presented in Section [3] followed by an analysis of the com¬ 
panion mass-ratio distribution in the overlap region of these 
two populations and a discussion of the BD desert and its 
claimed consistency with empirically determined IMF mod¬ 
els ( |Reggiani & Meyer|201 1 in Section [4] In addition, the 
term peripheral fragmentation is introduced as the main for¬ 
mation channel of BDs. 


2. METHOD 

The basic idea is to quantify the deficit of analytically 
derived mass functions with respect to observationally con¬ 
strained IMFs. The resulting deficit defines the RMF. 


2.1. Observational IMF models 


The stellar IMF is based on an extensive analysis of obser¬ 
vational data from young stellar clusters and in the Galactic 
field ( |Kroupa et al.|1993]|Kroupa|2001| i. Because the major¬ 
ity of multiple systems remain unresolved such IMFs need to 
be interpreted as system mass functions unless a careful nu¬ 
merical binary correction is used. 

The empirical IMF by C05 used here is a system IMF based 
on the earlier IMF by |Chabrier] ( j2003) >, assuming BDs and 
stars forming the same way, thereby simply occupying differ¬ 
ent regions in the same continuous mass function, an assump¬ 
tion also made in the past by|Kroupa|i 
l indr 


(J200T. 2002) TK07, on 
body IMF that has been 


the other hand, provide an individual-body 
transformed into the corresponding system IMF by Monte 
Carlo random pairing among the star-like and the separate 
BD-like population introduced in TK07. This accounts for 
the observational evidence for two separate (albeit related) 
formation channels. In particular, there is a lack of observed 


BD companions to stars, especially for small separations (Me 
Carthy~et al.|2003[ Grether & Lineweaver|2006 1 , whereas the 


statistical properties of binary BDs differ largely from stel- 
lar ones (|Bouy et akl2003 Burgasser et al.|2003l|Martm et al.| 


|2003||Close et al.|2003). As argued in|Thies & Kroupa|(|2007|!, 

the different binary properties of BDs and VLMSs on the one 


hand, and stars on the other, consequently suggest the exis¬ 
tence of two separate formation channels. This thus leads to 
the requirement for two separate IMFs, each corrected for un¬ 
resolved binaries, to describe the real individual-body mass 
function of a star-forming event. 

These IMF descriptions, however, are purely observation- 
ally motivated and therefore do not explain the underlying 
physical processes leading to the actual mass spectrum. 


2.2. Analytical IMF models 

In recent years, several attempts have been made to under¬ 
stand the physics behind star formation and to reproduce its 
outcome. This section deals with the analytical models by 
PN02 and HC08. In both PN02 and HC08 the prestellar clump 
mass function is described by an analytical function with the 
Jeans mass, the length scale of the initial inhomogeneities and 
the Mach number, Ml, as the main parameters (see Equation 
(24) in PN02 and Equation (44) in HC08). Both studies de¬ 
rive the IMF from the mass distribution of gravitationally un¬ 
stable clumps based on empirical power spectra of turbulent 
flows within the molecular clouds. While the particular for¬ 
malism differs in PN02 and HC08, fragmentation due to the 
supersonic interaction of gas sheets is the engine for forming 
Jeans-unstable clumps in both models. The resulting prestel¬ 
lar clump mass function is then transformed into a stellar IMF 


assuming a star-formation efficiency of 30 %-50 %. The mod¬ 
els have been tested in the range 3 < Ml < 20. In this paper, 
we restrict the model to the Mach number AJ = 12 in the case 
of HC08 (Figure 5 top therein) and Ml = 10 in PN02 because 
these Mach numbers are suggested as the most likely values in 
these papers (see p. 873 in PN02 and p. 406 in HC08). With 
these values both models match each other well in the stel¬ 
lar regime despite their different formulations. Higher Mach 
numbers would mainly shift the mass spectrum toward lower 
masses and decrease the MF toward the high-mass end. As 
shown in Figure [T] the PN02 Ml = 10 mass function (long- 
dashed curve) and the Ml = 12 HC08 mass function (solid 
curve) closely match in the stellar regime but are slightly dif¬ 
ferent at the low-mass end. The resulting mass function de¬ 
viates significantly from the observed stellar+substellar mass 
function in the substellar mass regime, as shown in Figures 
[2] and [3] for both theoretical models compared to the observa¬ 
tionally constrained IMFs by C05 and TK07. To quantify this 
difference the model IMF, either HC08 or PN02, is subtracted 
from an observational reference mass function that is taken 
from C05 and TK07 in order to determine two estimates of 
the RMF: £ res , 


£res(tti) — CobsC^t) £theo(^l) , (1) 

where G hco is any of the theoretical IMFs and £ 0 b s refers to 
any of the observationally constrained IMFs. In general, a 
mass function £ is defined as the differential number N over 
the differential object mass m: 

AN 

£('«>=—, ( 2 ) 

Am 

and, in the logarithmic scale 


£l ( m ) ■■ 


AN 


d log 10 zh/M© 


TlnlO )m/M©£(m). (3) 


The TK07 IMF is the canonical IMF (Kroupa et al. 2013 i 
that takes into account that BDs and some VLMSs" need to 
be added as an additional population called BD-like, whereas 
most stars belong to the star-like population: 


£ B D('»)-^p 0 pk 
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where 7?. pop is the 

BD-like to 

star-like population 

ratio 


0.5 


-1.3 


(T^pop « 0.2, TK07) and k m = V007 , 
of the stellar IMF at 0.5 M. . Here, m, 
m. 


■Mgci relation (Weidner et al.|20f0| |2013| Equation 
Pflamm-Altenburg et al. 2007) and approaches 150 M. 


ensures continuity 
follows from the 
( 10 ) 
z© Ra¬ 
the most massive clusters. Note the overlapping mass ranges 
of the populations indicating that bodies between about 0.07 
and 0.15 M© may belong either to the star-like or the BD-like 
population. 

At the high-mass end of the BD-like population, the sharp 
truncation used in TK07 has been replaced in this study by 
a steep power-law function to reduce numerical artifacts in 
the sum IMF. We chose a power-law exponent of 10 to keep 
the effect on the BD-like to star-like ratio negligibly small. 
There is no such declining function applied to the star-like 
population that is intrinsically smoother due to the dyna mical 
population synthesis (DPS) method described in Section [23] 

































3 



l°9i 0 m/M e 


Fig. 1.— Comparison of the analytical system IMF by PN02 for A4 = 10 
(their preferred value, dashed curve) with that by F1C0B f.Vl = 12, solid 
curve). Both functions are scaled in this plot for equal peak height for com¬ 
parison only. Although they nearly match in the stellar regime there are slight 
deviations in the substellar regime as well as in the positions (i.e. masses) of 
the peak. 


In addition to this, star-like objects below 0.06 Mg and BD- 
like ones above 0.3 Mg are not considered. It should be noted 
here that the theoretical mass function obtained by Thies et al. 
( j2010[ ) from smoothed particle hydrodynamics (SPH) simula- 
tions also shows a steep decline above 0.1 Mg rather than an 
exact truncation. 

Since the C05 stellar IMF used by HC08 is a system mass 
function rather than an individual-body-mass function, the 
TK07 canonical IMF has also been transformed into the corre¬ 
sponding TK07 system ma ss fu nction using the Monte-Carlo 
model described in Section [273] The stellar component of the 
canonical TK07 system IMF is normalized to the C05 system 
IMF in the whole mass regime (0.01-150 Mg), i.e. 

150 Mg 150 Mg 

J ^TK07,sys(w)dm= J £c05,sysO»)d'M • (5) 

0.01 Mg 0.01 Mg 

Here it is assumed that in the C05 IMF stars and BDs share 
the same IMF without a discontinuity or overlap, so the mass 
regime above 0.08 M 0 corresponds to the stellar component 
in the TK07 model. 


The BD multi 

tlicity fraction is only about 10%-20% 

(Bouy et al. 

2003, 

Close et al. 2003 Martin et al. 2003 Kraus, 

White, & H 

illenbrand 2006 

and Law et al.1120081 and is as- 


sumed to be equal for both C05 and TK07 IMFs. 


Further, it has to be noted here that the RMF is based on 
the mass range between 0.01 and about 0.3 Mg (log 10 m/Mg 
between -2 and -0.5) only, because the residuals in the stel¬ 
lar regime are due to slightly different functional descriptions 
rather than being physically constrained and are therefore ne¬ 
glected here. The residuals between the observational stellar 
IMF fits of C05 and TK07 and the HC08 and PN02 models 
are of potential interest to later work, but they are not further 
considered in this paper. 

2.3. Monte-Carlo model of the binarity of stars and brown 

dwarfs 

Besides the mass function itself, the binarity as a function 
of the primary-star mass and the CMRD are also important 
characteristics of stellar populations. They are studied here 
using a Monte Carlo approach with star-like and BD-like ob¬ 
jects drawn from separate IMFs. 

The binarity or binary fraction, /, is defined as the ratio of 
the number of binary or higher-order systems, A4,j n , to the total 


number of systems, N sys . Here, the term system includes mul¬ 
tiple systems and singles (their number being noted as (V sng ) 
as well. Then 

^ bin (6) 






sys 


(V sng + A)-, M 


For the star-like population we apply the binary DPS method 
developed by Marks, Kroupa, & Oh ( |2011} and Marks &] 


Kroupa| (|201 1 1 , hereafter referred to as dynamical or DPS 
pairing. In DPS the binary stars are formed in a popula¬ 
tion of embedded clusters within which they are dynami¬ 
cally processed to yield the Galactic disk stellar single-plus- 
binary population. An attractive feature of this theory is its 
underlying assumption of the universality of binary proper¬ 
ties of late-type stars being consistent with observational data 
(Marks & Kroupa 2012 Leigh et al. 2014| l. For initial bi- 
naries with intermediate to large separations the DPS pair¬ 
ing method applies random pairinqjbelow a primary mass of 
5 Mg and ordered pairing (such that q > 0.9) above. Here, 
q = m comp/'Mpnm < 1, where m comp is the companion mass 
and Wprim is the mass of the primary star. This initial bi¬ 
nary population is then altered by dynamical evolution. Close 
binaries with orbital periods below about 10 days undergo 
eigenevolution ( |Kroupa|1995| > and tend to equalize the com¬ 
panion masses. Note that this eigenevolution term alters the 
very-low-mass end of the starlike IMF. For the purpose of this 
work, however these effects only play a negligible role. Here, 
the initial or primordial binary fraction is 100%, i.e. it is 
assumed that all stars form in binaries. The final (after dy¬ 
namical processing in the embedded cluster) overall binary 
fraction is about 40 % (i.e. / = 0.4) but varies as a function 
of the primary-star mass. For M-dwarfs, in particular, it is as 
low as 25 %while G-dwarfs show about 56 % binarity. The 
binary fraction approaches 90% for O stars. For the BD- 
like population we chose an overall binary fraction of 20 % 
(i.e. / = 0.20), in accordance with TK07, TK08. About half 
of the members of observed average stellar populations are 
binaries, most of them remaining unresolved in typical star- 
cluster surveys. However, very young and likely dynamically 
unevolved populations like the Taurus-Auriga association ex- 
hibit almost 100 % binarity ( Kroupa et al.|2013) |Duchene~&] 


|Kraus||2013| |Reipurth et al.| 2014| l. The number of systems 

must not be confused with the number of individual bodies, 
/V h(K i = (V sng +2(Vbi n . Since higher-order multiples are relatively 
rare ( Goodwin & Kroupa|2005f they are summarized within 
the binary population in this work, so the total number of bod¬ 


ies is 


(Vbod — —(Vbin A SI 


(7) 


The CMRD describes the relative number of binaries as a 
function of the companion-to-primary mass ratio. Observa¬ 
tions reveal a continuous decline of / as a function of the 
primary-object mass which has been interpreted as a contin¬ 
uous transition from the stellar to the substellar re gime flJo-l 


ergens 2008] |Kraus & Hi 1 lenbrand| |20 1 2[ but see Thies & 
Krou pa|2008[ l. There is also a shift to wards more eq ual-mass 


binaries (q = 1) for VLMS and BDs (Dieterich et alT| |2012| >. 
These properties of the stellar population are well reproduced 
by DPS such that the origin and properties of binary popula¬ 
tions are well understood. 


- Marks et al. (2014, in prep.) show that random pairing with subsequent 
dynamical processing does indeed repro duce at least the observed low-mass 
stellar population (see also|Kroupa| 1995 j. 
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In this study we used a Monte Carlo approach with the 
TK07 IMF (Equation ([4])), i.e., with two separate populations, 
BD-like and star-like, ranging from 0.01 to 0.15 Mg and from 
0.06 to the maximum stellar mass of 150 Mg, respectively. 
The slightly larger mass ranges compared to Equation 0 are 
due to the steep power-law decline added to the mass bor¬ 
ders to reduce numerical artifacts. The masses are drawn ran¬ 
domly from each IMF where the relative normalization be¬ 
tween both populations is simply obtained by the number of 
objects drawn from each IMF. 

BD-like objects are assumed to form as single substel- 
lar cores some of which are subsequently paired to binaries 
within a dense dynamically preprocessed environment like 
ive extended accretion disk (S t am ate 11 os et al.|2007 a, 


a massive 


Thies et al. 2010 Basu & Vorobyov||2012| i. Similarly, stel¬ 

lar binaries are assumed to be assembled from individually 
formed and subsequently paired stars. After the individual- 
body populations have been drawn from the IMF, two differ¬ 
ent methods are used for assembling the binaries. For stars 
we apply the DPS pairing method mentioned above. 

BD binaries, on the other hand, are created by first draw¬ 
ing two objects from their separate population characterized 
by £bd (Equation |4j) with the pairing probability /? pair being 
determined by a power law. 


Impair “ COnSt tj , 


( 8 ) 


following the power-law bias scheme used by |Goodwin| 
( |2013| ) for the mass ratio distribution of binaries from sec¬ 
ond (i.e. binary-forming) fragmentation. Whereas the case 
7 = 0 corresponds to random pairing from the IMF, 7 > 0 de¬ 
scribes a biased pairing rule with an increasing pairing prob¬ 
ability toward equal-mass components. The extreme case 
of perfect equal-mass pairing corresponds to 7 = 00 but is 
in disagreement with data from the Very-Low-Mass Binary 
Archive (VLMBA; Burgasser et al. |2007| which also con¬ 
tains a few unequal substellar binaries. In practice, first an 
array of BD-likes is generated by drawing randomly from £bd 
(E quation [4]). All subsequent pairings to make binaries are 
performed on the array only. Random pairing is performed by 
randomly drawing a companion from the array for each ob¬ 
ject of the same population. If the companion is more massive 
than the considered object it becomes the primary, and other¬ 
wise it is the secondary component. In the case of biased pair¬ 
ing an additional decision is made whether a companion ran¬ 
domly drawn from the population is accepted to be a partner 
or rejected, depending on the probability p plllr in Equation (fsl. 
A rejected partner may later be chosen as a partner to another 
object. Objects that, on the other hand, are already bound in 
a binary are skipped in the pairing procedure henceforth. The 
biased pairing procedure is iterated until the required binary 
fraction / is achieved. Because the objects are chosen in ran¬ 
dom order, this does not introduce any additional bias besides 
Ppair t° the binary populations. The method ensures that the 
slope of the canonical IMF is retained, in contrast to methods 
that select components according to a mass-ratio distribution 
only. 


3. RESULTS 

3.1. Residual Mass function for Semianalytical Models 

In this section we present the RMFs obtained from the 
BD/VLMS deficit of the analytical IMF models by PN02 and 
HC08. The RMF is derived for a particular and, according to 
PN02 and HC08, typical parameter set, in particular the tur¬ 


bulent Mach number. Other parameters yield different RMFs, 
so only the general functional shape for a likely parameter 
set is presented here. It is further shown that such an RMF 
provides a correction term to the stellar IMFs in order to fit a 
reasonable overall IMF without the BD deficit observed in the 
purely stellar ones. Figure [2] shows the results of our calcu¬ 
lations for the analytical model from HC08 with the C05 and 
TK07 system IMF as the observational reference (upper and 
lower panels, respectively). In both cases the RMF is lim¬ 
ited to the mass range below 0.3 Mg (log 10 w/Mg < -0.5) 
because the slightly different shapes of the IMFs in the stellar 
region are beyond the scope of this paper. The local min¬ 
ima near log 10 w/Mg = -0.8 (upper panel) and -1.4 (lower 
panel) occur because the difference between the HC08 mass 
function and the empirical mass function almost vanishes lo¬ 
cally. This behavior is highly sensitive to the normalization 
of the IMFs and is most prominent in the case of the discon¬ 
tinuous TK07 IMF. Peaking near the hydrogen-burning mass 
limit the RMF declines steeply and effectively vanishes above 
log 10 ih/Mq = -0.5. The peak mass, however, is better repre¬ 
sented in the residual to TK07 in Figure [2] whereas the RMF 
from C05 appears to be shifted toward lower masses by about 
a factor of two (i.e. an offset of about -0.3 on the logarithmic 
scale). 

In the case of the PN02 analytical model, as shown in Fig¬ 
ure [3] the result with respect to the C05 IMF covers a mass 
range similar to HC08 but is also continuous. This is ex¬ 
pected because both analytical models have a functional form 
similar to the C05 IMF, namely an extended log-normal-type 
shape. In the case of the TK07 system IMF (Figure [3]) there 
is a prominent “dip” in the RMF around log Io m/M 0 = -1.3 
(i.e. m = 0.05 Mg). The reason is the near-equality of TK07 
and PN02 at this point, which varies sensitively with the nor¬ 
malization of the TK07 system IMF and its substellar compo¬ 
nent. As with the HC08 model, the RMF of PN02 versus C05 
peaks at lower masses by a factor of about two. 


3.2. Comparison with Simulation Data 

The RMF obtained in the previous section is compared here 
with the mass spectrum of substellar clumps formed in the 
SPH models of Thies et al.~j (j2010| l. Therein the induced frag¬ 
mentation of massive extended circumstellar disks due to per¬ 
turbation by passing stars in an embedded cluster has been 
studied. Using an SPH code with a radiative cooling approx¬ 
imation (Stamatellos et al. |200 7b ) gravitational instabilities 
have been demonstrated to form through tidal perturbations 
and to form compact clumps with typical masses between 
0.01 and 0.15 Mg. The mass function of the thus-formed 
clumps, the “SPH mass function” (SPH MF) hereafter, has 
been shown in Thies et al. (1201 Oil to be in agreement with 
the substellar component of the observed system IMF from 
TK07. In the current study, the SPH mass function is based 
on 80 objects formed in 29 computations, some of them per¬ 
formed in continuation of TK07. 

In Figure [4] the (scaled) SPH mass function is compared to 
the corresponding RMF. The overall shape of the RMF essen¬ 
tially matches the SPH mass function which features a gen¬ 
eral increase with increasing mass, a peak around 0.1 Mg and 
a rapid decay for higher masses, thus characterizing a pop¬ 
ulation that is restricted to the substellar and very-low-mass 
stellar regime. As mentioned by Thie s et al.|(|2010 1 , these re¬ 
sults a re quite similar to those in Stamatellos & Whitworth 
(2009j ) who computed BD formation in self-fragmenting 
disks. Therefore, the results from triggered fragmentation can 
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FIG. 2.— Upper panel: the analytical IMF model for Ad = 12 by HC08 
(solid line) compared to the empirical IMF by|Chabrier|{2005) (dashed line). 
These functions, originally defined as system mass functions, have been nor¬ 
malized to match in the stellar regime. The dotted line represents the residual 
mass function, i.e. the difference between both mass functions. The gaps in 
the RMFs are caused by the intersection of the HC08 model IMF with the 
empirical IMFs. They are highly sensitive to the normalization of the IMFs. 
Lower panel: sa me as in the upper panel but with the combined bimodal 
IMF according to|Thies & Kroupa|{2007| (dashed line). Note that the RMF 
is truncated for log > -U.5 because the different functional shapes 

of the stellar parts are not considered here. 



Fig. 3. — Upper panel: the same as in Figure[2]but with the clump mass 
function model by PN02 instead of HC08. Lower panel: same as Figure[2l 
lower panel, but with the clump mass function model by PN02 instead of 
HC08. Here, the RMF does not contain any gaps. 


be taken as representative for the general disk fragmentation 
outcome. Both scenarios assume that the substellar clumps 
formed in a fragmenting disk are ejected by either mutual dy¬ 
namical interactions of clumps in the dis k or in subsequent 
stellar encounters. According to |Basu & Vorobyov| (|2012l 
some of these clumps are disrupted by the tidal forces during 
the ejection or are destroyed because of migration onto the 



-2 -1.5 -1 -0.5 0 

log io ™/M Q 

Fig. 4.— Residual mass functions of HC08 (upper panel, from FigureUl 
and PN02 (lower panel) with respect to the system mass functions by C05 
(dashed line, from Figure [5} and TK07 (dotted line) as well as the aver age 
of both (solid line) in comparison to the SPH mass function from|Thies et al.| 
(20TO] including post-population SPH data), shown as a histogram (normal- 
lzecTin log 10 £z, to the RMF). 


central star. The impact of this effect on the mass distribution 
will be subject to future research. 

Given the uncertainties of the PN02 and HC08 MFs the 
RMF based on them can be considered to be in agreement 
with the substellar IMF in both TK07 and |Thies et al.| ( [20T0| . 
Interestingly, the average of both RMFs (each weighted as 
50 %) provides an even better fit to the SPH data (see the solid 
curves in Figure[4]). Use of such an average may be motivated 
by the expectation that the true IMF below « 0.2M 0 is likely 
between the TK07 and C05 quantifications. 


3.3. Monte Carlo Study Applied to Observations 

The IMF from the Monte Carlo study described in Section 
|2.3| is shown in Figure [5] in comparison to the system mass 
function (dash-dotted curve), the primary-body mass function 
(dashed curve) and the individual-body mass function (solid 
curve). The upper panel displays the sum IMF with both pop¬ 
ulations combined, as it would be seen by an observer. In the 
lower panel, both populations are plotted separately. 

We note that Figure 12 in Kirk & Myers | ( |2012| >, which 
shows the observed average mass functions of young stel¬ 
lar clusters, groups and isolated objects across the stellar- 
substellar boundary, reveals a prominent step near the stellar- 
substellar border that is very similar to the one that emerges 
from our analysis as plotted in the upper panel of Figure [5] 
The system MF (dash-dotted curve) is to be compared with 
the mass functio n by |Kirk & Myers | (|2012|). Because |Kirk| 


& Myers| ([201 2) only studied low-densrty, dynamically un¬ 
evolved populations, the stellar system IMF is modelled here 


with a binary fraction of 80 % (see also fig. 10 in [Marks et~ak 


|2011| . For the Galactic field, for which the rest of our study 
has been performed, the binary fraction is only about 40-50 % 
for low-mass stars. 

Figure [6] depicts as the solid line the binarity as a function 
of the primary-body mass (i.e. the system mass in the case of 
single objects) from the DPS pairing model for the Galactic 
field population for embedded clusters with a half-mass radius 
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Fig. 5.— Upper panel: the combined IMF of BD-like and star-like objects 
from our Monte Carlo calculations for individual-body masses (solid curve), 
primary-body masses (dashed curve), and system masses (dash-dotted curve). 
The BD-like binarity is 20 % here, and the star-like one is here about 80 % on 
average, to match the case of dynamically unevolved populations (se e tex t). 
For comparison, the histogram for grouped stars in |Kirk & Myers| {2012| 
"KM12”), which also shows a discontinuity near the stellar-substellar border, 
is included as the dotted histogram. Lower panel: the IMFs for the BD-like 
and star-like populations plotted separately. The line patterns are the same as 
in the upper panel for each separate IMF. For comparison, the individual-body 
sum IMF is superimposed (thin dotted curve, the same as the solid curve in 
the top panel). Note that the very-lo w-mas s end of the star-like IMF is slightly 
depleted due to eigenevolution {Marks et al.|201 1). 


of 0.1 pc for the underlying cluster population. The long- 
dashed line indicates the case of a half-mass radius of 0.2 pc 
which yields an up to 20 % higher multiplicity. These mode l 
Galactic field populations, computed with DPS (Section [23} , 
assume that all stars form in binary systems in a population 
of embedded star clusters that evolve and disperse their stars 
into the field. 

These two cases are not to be taken as statistically strict 
confidence limits but rather as two l ikely cases that have been 
studied in [Marks & Kroupa| (| 20 TT), and the initial stellar bi¬ 
nary population is indicated by the shaded region at the top. 

The stellar binary fraction resulting from DPS is about 
/ = 0.4 on average in the low-mass region below 1 M 0 for 
the case of 0.1 pc. The single-hatched region indicates the 
large uncertainty of the observed binarity in the substellar re¬ 
gion. For comparison. Galactic-field observational data are 
overplotted for A stars (VAST survey, De Rosa et al. 2014 as¬ 
terisks), G and late-F stars ( Duquennoy & Mayor|1991| filled 
squares), M to G stars (|Kroupa et aE 2003| open squares), 
M dwa rfs (|Fischer & Marcy||1992| upright triangle), and L 
dwarfs ( |Reid et al.|2008 upside-down triangle). 

Despite consisting of two separate populations (BD-like 
and star-like) the binary fraction declines continuously from 
the stellar toward the substella r regime. As d iscussed in 
Marks & Kroupa| ( |2011j l and Marks & Kroupaj ( |2012} , this 
is primarily due to the dynamical processing of the star-like 
population. In part, the mass overlap of the star-like and the 
BD-like population also contributes to this transition. It re¬ 
mains to be shown in future work that the apparent observed 
trend toward narrow binary separation distributions for late- 
type M dwarfs reported by Janson et al. ( 2014} can also be 
reproduced with DPS by this overlap. 


3.4. Binarity and Companion Mass-ratio Distribution 

The observable BD CMRD, which is the CMRD for all pri¬ 
maries (BD-like and star-like) below 0.08 Mq, is shown in 
Figure [7] As can be seen the case of strongly biased pair¬ 
ing with 7 = 2 (solid line) matches the measured data from 
the VLMBA better (and also see |Dieterich et al.||2012| ) than 
the random-pairing q { 7 = 0 , dotted line) and weakly biased 
pairing models ( 7 = 1 , dashed line). In addition, the case 
7 = 4 ( |Duchene & Kraus||2013[ dash-dotted line) is shown. 
The binarity of BDs near the hydrogen-burning mass limit of 
0.08 Mq is quite low in the model, only about 11 %- 12 % 
which is largely due to the “dip” in the binarity in the VLMS 
region of the DPS pairing method (see Figure [ 6 }. This lo¬ 
cal minimum is a direct consequence of the two overlapping 
populations. |Reid et aI7[ ( [2008} indeed report a similar binary 
fraction of 12. 5+3 % for L dwarfs. If a local minimum is con¬ 
firmed by future surveys with sufficient precision in mass this 
would further support an overlap of two separate populations 
in this mass region. However, in the current data this feature 
is within the uncertainties of the observational data. 


3.5. Contribution of Peripheral Fragmentation to BDs 


One important implication of the model discussed in this 
paper is the hybrid nature of BDs and VLMSs. Although 
there is a minor contribution by the star-like population from 
direct fragmentation, the majority of BDs are contributed by 
the BD-like population through dynami cally preprocessed gas 
(e.g. fragmenting circumstellar disks (fS t a m ate 11 os & Whit-| 
worth [2009] Thies et al.||2010 i. Because this preprocessed 
material often occurs in the peripheral regions of star forma¬ 


tion (e.g. in the outer parts of disks) we propose the term pe¬ 
ripheral fragmentation for the BD-like formation channel. Of 
all BDs between 0.01 and 0.08 Mq 64 % are contributed by 
the BD-like population, and 19% of M dwarfs between 0.08 
and 0.45 Mq are BD-like. This result, however, is highly sen¬ 
sitive to the chosen lower mass border of the star-like regime. 
Here, we assumed it to be 0.06 Mq. If, on the other hand, 
a sharp truncation of the star-like regime at 0.08 Mq is cho¬ 
sen, the BD-like fraction of BDs is 100 % but still 18 % of M 
dwarfs are BD-like. Because star-like BDs and VLMSs may 
be detectable by larger circumsubstellar disks or, if applica¬ 
ble, wide binary separations, future high-resolution observa¬ 
tions will help to further constrain the BD-like and star-like 
mass borders. 


4. DISCUSSION 

In this paper we have introduced the RMF as a correc¬ 
tion term for the analytical star-formation models by PN02 
and HC08 to match the observationally constrained IMF by 
Chabrier (2005 i and Thies & Kroupa| ( |2007 1. The effec- 
tive deficit of BDs and VLMSs in these theoretical models 
with respect to the observed BD statistics suggests the re¬ 
quirement for additional formation channels in these analy¬ 
tical IMF models. HC08 speculate that this deficit, being at 
least at the edge of significance (see uncertainties shown in 
Figure 4-23 in Kroupa et al.|[2013} , may be solvable by an 
improved algorithm that accountsTbr the effects of turbulence 
and other dynamical processing of the prestellar gas. The ex¬ 
pected contribution of these additional effects to the IMF may 
be understood as a separate population or, at least, as an ad¬ 
ditional formation channel of BDs and VLMSs and can for¬ 
mally be described by an additional correction term. It is also 
applicable as a test for future star-formation models. 
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Fig. 6 .— Binary fraction for BDs and stars as a function of the primary-object mass. The solid line shows the result of our Monte Carlo computation for 
stars via DPS pairing for the field population based on the integration of the star-formation outcome over all cluster masses. A half-mass radius r/, = 0.1 pc is 
assumed for all embedded host clusters. The average field-star binary fraction is / 0.5 for low-mass to intermediate-mass stars, in accordance with TK07. The 

long-dashed line refers to the DPS model with = 0.2 pc, corresponding to a higher binary fraction. The short-dashed curve represents the BD-like population 
with biased pairing. The single-hatched region in the lower left comer indicates the uncertainty of the substell ar binary fraction in the observational data. The 
smoothly shaded area at the top refers to the initial stellar binary population with a binary fraction of near 100 % (Marks & Kroupa|20lT) . Observational field-star 
data are overplotted for A stars (VAS T survey, De Rosa et al. 2014] asterisks), G and late-F stars (Duquennoy & Mayor 1 991| tilled squares), M to G stars fjKroupa 
et al. 2003 open squares), M dwarfs (Fischer & Marcy 1992 upright triangle), and L dwarfs (Reid et al. 2008, upside-down triangle). 




0 0.2 0.4 0.6 0.8 1 


Q - ^comp^prim 

Fig. 7.— Number distribution f(q ) for binary BDs between 0.03 and 
0.08 Mq as a function of the companion-to-primary mass ratio. The solid 
line shows the result of our computation for 7 = 2 , and the dashed and dotted 
lines refer to the case 7=1 and random pai ring, respectively. For co mplete- 
ness, also the more extreme case of 7 = 4 jDuchene & Kraus|2013) is also 
shown as the dash-dotted line. The histog ram shows the observa tional data 
from the Very-Low-Mass Binary Archive (Burgasser et al.|2007). The peak 
at q = 1 is primarily due to the contribution oFbinaries from the low-mass end 
of the star-like population because any binary with both components near the 
lower mass border of a population can only have nearly equal component 
masses. The area below each curve and below the histogram is normalized to 
1. 


This correction term has been identified and quantified here 
as the RMF which has a similar shape for both analyzed obser¬ 
vational reference IMFs by C05 and TK07. The most striking 
outcome is a general agreement between the RMF (Equation 
0 ) as obtained from the theoretical IMF and the observa¬ 
tional IMFs, and the mass function of fragments from SPH 
simulations ( |Thies et ak1|2010| l. In other words, both PN02 
and HC08 IMF models appear to describe essentially the pop¬ 
ulation of stars without a significant fraction of BDs. They can 
be considered to be successful models of the direct fragmen- 



log 10 m/M e 

Fig. 8 .— Comparison of the |Thies & Kroupa| (2007), |Chabrier| (2005a. 
and the extrapolated Bochanski et al.| , 21 ) I (I, IMF, the latter being used by 
Reggiani & Meyer ,20 1 Ithe original Bochanski IMF has been derived 
from data between 0.1 and 0 .8 Mq . Note the steeper decline in the Bochanski 
IMF below its valid mass range relative to the Chabrier IMF. 


tation process in molecular clouds. These models cannot, by 
their nature, accommodate the peripheral fragmentation, e.g., 
in accretion disks aroun d pr otostars, which yields the BD-like 
population (see Section [XT]). The most natural way to explain 
these results is, to the best of our knowledge, a composite 
population consisting of a star-like and a BD-like component 
such that a significant fraction of the BD part comes from an 


additional process acting during star formation (Reipurth & 

Clarke 2001, 

Stamatellos et al. 2007a, Thies et al. 2010, Basu 

& Vorobyov 

2012). Essentially this result follows for BDs 


being lesslikely to form directly in a star-forming molecular 
cloud than from dynamically preprocessed material. This is 
because in order to form only a BD a cloud core needs to have 
a high pre-collapse density (for it to be gravitationally unsta- 
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Fig. 9.— Snapshot from SPH model X002 by|Thies et al.|(2010|l showing 
substellar clumps (with sink particles marked with white points) around a 
fragmented accretion disk, all of them surrounded by their own accretion 
envelopes. Note that the escaping object (0.013 Mg, i.e. a low-mass BD) 
to the right has retained its accretion disk (visible as diffuse structure around 
the sink particle) even after dynamical ejection. The half-mass radius of the 
envelope is about 10 AU, and its mass is 2- 10 4 Mg, i.e. almost 2% of the 
mas of the escaping object. 


ble) while having no further supply of gas in order to limit 
its further growth in mass. This has already been emphasized 
by |Adams~& Fatuzzo| (1996| ). While there has been a recent 
discovery of a possible proto-BD in Oph B-l 1 with a mass of 
0.03 Mq (Andre et al.|2012| its further evolution and possible 
additional accretion from the surrounding molecular gas re¬ 
mains unclear. And, the existence of the BD-like population 
formed through peripheral fragmentation does not exclude the 
formation of BDs via the direct fragmentation channel, i.e. as 
star-like objects as is evident from the analytical PN02 and 
HC08 IMFs. 


4.1. Is There a Separate Substellar Population? 

There is an ongoing discussion whether a separate popu¬ 
lation is really needed to explain the observed properties. A 
similar discontinuity between planets and BDs is widely ac¬ 
cepted. |Chabrier et al.| (|2014|l even report a possible mass 
overlap of the BD and the giant planet regime as well and 
suggest a distinction between BDs and planets by their for¬ 
mation history rather than by deuterium burning. This can 
be seen as an analog to the separate substellar population ad¬ 
dressed by our work. jJumper & Fisher|(j2013|l claim to be able 
to model the trend toward more tightly bound binaries in the 
BD regime. However, because this study already assumes the 
masses to be drawn from a continuous IMF it a priori excludes 
the option of a two-component IMF. The AstraLux survey re¬ 
ports a narrower binary separation distribution for M-dwarf 
binaries compared to binaries with FGK primaries (|Janson| 
et al.|2012 and Janson et al.|2014 respectively). Such a nar- 
row peak in the separation distribution may be related to sen¬ 
sitivity artifacts as well as, at least for late-type M dwarfs, due 
to the mass overlap with the BD-like population that extends 
up to about 0.2 Mq. _ 

In contrast, in a study of dynamical binary evolution Parker 
& Goodwin| (]2011J) found that the field populations of M 
dwarfs and very-low-mass binaries must reflect very different 
birth populations, since their dynamical processing is essen¬ 
tially the same. 


4.2. The Brown Dwarf Desert 

Similarly, the BD desert itself has been disputed by several 
groups who suggest that it is more related to the companion 
mass ratio rather than to the absolute mass. However, a recent 


survey by Dieterich et al. (2012) further supports the existence 
of a true BD desert, i.e. a dearth of substellar companions to 
primary stars ranging from M to G dwarfs. They found an 
effective lower mass limit of stellar companions near 0.1 M 0 
which cannot be related to incompleteness but rather points 
toward a discontinuity in the pairing statistics close to the 
hydrogen-burning mass limit. They particularly point out that 
this discontinuity is essentially independent of the primary 
star mass. This clearly supports a separate substellar popula¬ 
tion contributing the majority of BDs and is also in agreement 
with SPH studies of whole star-forming clouds by, e.g.,|Bate| 
et ah] (2003]) and (Bate] (2009a) who reproduce the formation 
of BDslargely from fragmenting circumstellar disks. They, 
however, point out that these simulations tend to overproduce 
BDs and attri bute this to the lack of radiative heating in their 
model (Bate|2009b| |201 2) , It should also be noted that [Bate] 
(2012) could reproduce a mass function in good agreement 
with the C05 IMF and that the previous overproduction of 
BDs was avoided there because the radiative heating reduced 
the frequency of disk fragmentation. 

In a related Monte Carlo study performed here, we also 
found a good agreement of the two-population composite 
model with observational data on field very-low-mass stel¬ 
lar and BD binaries (Figure [6]. The binarity is well repro¬ 
duced down to the stellar-substellar border as a continuous 
function of the primary-star mass in agreement with the ob¬ 
servations. More importantly, the mass-ratio distribution of 
very-low-mass binaries deduced from observational data is 
well reproduced in our model (Figure [7J, once more support¬ 
ing a bimodal star and BD formation scenario. Whereas star- 
like binaries are well represented by DPS (Marks et al.|2015(), 
the BD-like binaries apparently do fit a biased pairing where 
the probability of pairing rises toward more equal component 
masses following Equation (8]). Here, the cases 7 = 0 (sim¬ 
ple random pairing) and 1 < 7 < 2 (biased toward equal-mass 
binaries) were compared. The biased pairing case is in bet¬ 
ter agreement with the observational data. It remains to be 
studied whether postformation dynamical processing through 
stellar-dynamical encounters in their birth embedded clusters 
is responsible for this apparent preference of more equal-mass 
binaries over more unequal ones in the BD mass range. How¬ 
ever, both simple and biased pairing are in reasonable agree¬ 
ment within the uncertainties of the data if a low binary frac¬ 
tion of 10% within the BD-like population is assumed. A 
binary fraction of less than 10% in the theoretical disk frag¬ 
mentation outcome and the very narrow range of semimajor 
axes have also been found by Thies et al. (2010) , as well as 
the survival of circumsubstellar accretion disks, as shown in 
Figure [9] The overall binary fraction for systems within 0.03 
and 0.08 M 0 , irrespective of the population to which they 
belong, is around 20 % for both simple and biased pairing 
(i.e. the value chosen for the BD-like population), if the local 
minimum in the binarity function near the stellar-substellar 
boundary (Figure [6) is assumed to be an artifact, whereas it is 
as low as 11 % otherwise. However, observations of VLMSs 
indicate a binary fraction as low as 12.5 % ( Reid et al.|2008[ ), 
so this “dip” may indeed be real and thus be a result of the 
overlap of two separate populations. 

4.3. Companion Mass Ratio Distribution 

There is an ongoing dis cussio n whether random pairing is 
applicable for stars. Reggiani & Meyer ( 201l) claim to rule 
out this model in favor of a universal CMRD. However, their 
results only cover a narrow range of validity of the underlying 
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model IMF (as specified by Bochanski et al.j |2010 see Fig¬ 
ure [8]) that only covers M, K and late-G stars, but no BDs. 
Consequently, the CMRD derived from it is quite limited, es¬ 
pecially for M dwarfs, for which no hypothesis test on the 
pairing rule can be made on this basis. Reggiani & Meyer 
(201 l) perform pairing experiments that actually exceed the 


range of validity of the underlying mass function by Bochan- 


ski et al. ( 

O' 

© 

Goodwin 

(2013) 


suggests that a random partition of protostel- 
fragments is in better agreement with observations 
of stellar binarity than an initial binary population drawn from 
the IMF with subsequent dynamical and eigenevolution. Be 
cause that study was based on the IMF by Chabrier (2003 


without a separate BD treatment there cannot by any satis 
factory agreement with the observed substellar CMRD. Even 
more importantly, the initial binary population needs to be 
modeled and tested against dynamically barely evolved stellar 
populations like Taurus-Auriga rather than against the Galac¬ 
tic field. The field population then results from an integration 
of all dynamically evolved clustered populations following 
the embedded cluster-mass function ( Marks & Kroupa|2011 
|Marks et al.|2015) . 

As a final remark, it is being argued that the existence of 
disks around young BDs implies that these form like stars sup¬ 
posedly supporting the continuous IMF scenario. Figure [9] is 
one of many examples where a young BD with an accretion 
disk formed from peripheral fragmentation in a circumstellar 
disk is nudged away to become a free-floating BD with a disk. 

5. CONCLUSION 

As the main conclusion from the first part of this study we 
emphasize the difficulties of theoretical star-formation mod¬ 
els to describe both stars and BDs by a single mechanism, 
namely from direct cloud fragmentation. It remains neces¬ 


sary to treat BDs separately which implies a separate albeit 
related formation channel for the maj ority of BDs. This al¬ 
ready follows from the work of Adams & Fatuzzo (1996). 
The theoretical evidence by Stamatellos & Whitworth (2009), 


|Thies et al.| (|2010J), and |Basu & Vorobyov|(|2012|) supports a 
separate population by fragmentation of extended young cir¬ 
cumstellar disks. If O stars are close by, the accretion enve¬ 
lope of VLM protostars may be photoe vapo rated, leaving an 
unfinished substellar embryo ( Kroupa & Bouvier||2003 1. In 
addition, the embryo-ejection model by Reipurth & Clarke 
(j2001 () gives an example of BD formation by ejection of un¬ 
finished stellar embryos out of multiple protostar systems. 
Because these mechanisms are not covered by the analytical 
cloud fragmentation models by PN02 and HC08 they do not 
contribute to the resulting theoretical clump mass function. 
The same is true for BD formation in dense gaseous filaments 
that form due to the gravitational pull of surrounding stars 
( Bonnell et al.||2008 ) rather than being described by random 
density fluctuations assumed by PN02 and C05. Therefore, 
an analytical model also covering BDs and VLMSs must in¬ 
clude such separate formation channels via dynamically pre- 
processed material (disks, unstable multiple stellar embryo 
systems, photoevaporation, tidally shaped gas filaments and 
so on). The normalization may be adjusted empirically or de¬ 
duced from physical relations based on typical disk-to-stellar 
mass relations and by the application of turbulence and the 
Jeans criterion within this preprocessed material. The devel¬ 
opment of such an analytical or semianalytical model should 
be the subject of future studies. The observational data thus 
very strongly suggest that the dominant fraction (Section [33) 
of BDs and VLMSs results from peripheral fragmentation that 
contributes a separate IMF from that of stars. 

This project has been funded by DFG project KR1635/27. 
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